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Abstract 

Experimental data on the shape of hadronic momentum spectra are compared with theoretical 
predictions in the context of calculations in the Modified Leading Log Approximation (MLLA), under 
the assumption of Local Parton Hadron Duality (LPHD). Considered are experimental measurements 
at e + e~ -colliders of £*, the position of the maximum in the distribution of £ p = \og(l/x p ), where 
x p = p/pbeam- The parameter £* is determined for various hadrons at various centre of mass energies. 
It is interesting to look at the dependence of £* on the hadron type. This is used to study the influence 
of the hadron type on the cut-off scale Qo in the parton shower development. The dependence of £* on 
the centre of mass energy is seen to be described adequately by perturbation theory. The approach is 
made quantitative by extracting a value of a s (mz) from an overall fit to the scaling behaviour of £* . 







1 Introduction 

During a few years of LEP running, a large amount 
of information was collected on identified hadron 
species in jets. The higher centre of mass en- 
ergy of LEP, compared to past e + e" accelerators, 
also makes it easier to separate the behaviour of 
hadrons with a high momentum, which are corre- 
lated strongly to the primary quark, from those with 
a low momentum, created mainly during fragmen- 
tation and hadronisation. 

Clear scaling violations have been observed in 
the shape of the charged particle x p -distribution as 
function of the centre of mass energy. The strong 
coupling constant has been extracted from these 
scaling violations, using the behaviour at high mo- 
menta: 0.2 < x„ < 0.7. Due to the larger statistical 
errors this is not possible for individually identified 
hadron species. 

However, low momentum data for specific types 
of hadrons maybe be used to study the properties of 
jet-evolution and hadron-formation in the context 
of the LPHD hypothesis and MLLA calculations 
of parton spectra [1-4]. The assumption of 'Local 
Parton Hadron Duality' (LPHD) states that a cal- 
culated spectrum for 'partons' in a 'parton shower' 
can be related to the spectrum of real hadrons by 
simple normalisation constants. These constants 
have to be determined by experiment. A second 
assumption is that the low momentum part of the 
spectrum is not influenced in a significant way by 
hadrons that are correlated to the primary quark. 

Calculations of the parton spectra in the 'Mod- 
ified Leading Log Approximation' (MLLA) take 
into account next-to-leading logarithms in a con- 
sistent fashion. The physical mechanism relevant 
to these next-to-leading logarithms is the coherent 
emission of soft gluons inside a jet, leading to an 
angular ordering and an effective transverse mo- 
mentum cutoff for the partons. Parton jets develop 
through repeated parton splittings, resulting in an 
increase of the multiplicity at lower momenta. The 
interplay of coherent emission of gluons and the cre- 
ation of hadrons causes this spectrum to be cut off at 
very low momenta. Calculations predict the shape 
of the ^ p distribution. The resulting 'hump-backed' 
distribution is nearly gaussian. As an example, the 
Monte Carlo spectrum of the A baryon in Z° decays 
can be seen in figures 1 and 2 . In the following 
will be determined for various types of hadrons and 
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Figure 1: The dependence of the A spectrum on the 
species of its 'parent' particle in JET SET. Particles 
produced directly in the fragmentation process have 
the Lund-model string as 'parent'. The scales of 
both x p and ^ p are given for comparison. 
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Figure 2: The dependence of the A spectrum on the 
flavour of the primary quarks in Z° decays at LEP, 
as predicted by JET SET. Notice that the high and 
low momentum regions are not separated for A 's in 
s, c and b events. 
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Table 1: The values of ^* for various hadrons as determined from the momentum spectra as measured 
by e + e~ experiments at different values of \fs. For the i] meson, the L3 experiment [10] reported a value 
£ = 2.60 ±0.15. 



at various centre of mass energies. Subsequently 
a comparison is made with the theoretical calcula- 
tions. 

Figure 1 shows that LPHD is not an obvious as- 
sumption: in the JETSET Monte Carlo, the spec- 
trum of e.g. a A baryon depends on its 'parent' 
particle. The LPHD assumption states that the sum 
of these spectra is proportional to the spectrum as 
calculated for a parton shower with a correctly cho- 
sen energy cut-off. 

Moreover, the JETSET momentum spectrum of 
the A baryon (figure 2) shows a strong dependence 
on the flavour of the primary quarks. For the heavier 
quarks, s, c and b, the high momentum contribution 
correlated to the primary quarks is not separated 
from the low momentum contribution related to the 
fragmentation process. This can be expected to 
cause shifts in that depend on the jet-flavour. 

The spectrum of all charged particles at LEP has 
been fitted to the MLLA distribution, with free nor- 
malisation factors for pions, kaons and protons, de- 
pending on the centre of mass energy [5,6]. It was 
quite surprising [3] that the MLLA functions can 
also fit the high momentum part in the data for 
pions. This is probably a coincidence, since the 
heavier K° meson has a spectrum that can not be 
fitted as nicely at large x p . 

One way to see whether the observed shape of the 
spectrum is really due to coherent gluon emission 
is to compare experimental data to predictions of 



the JETSET Monte Carlo program, with the coher- 
ence (angular ordering) either turned on or turned 
off. The large dependence on the primary quark 
flavours and the fact that coherence in JETSET is 
not really necessary to fit the experimental data has 
led some to the conclusion [7] that coherence can 
not be demonstrated on the basis of hadron mo- 
mentum spectra. Here an attempt will be made to 
infer more information using the properties of iden- 
tified hadrons, concentrating on the observable £*. 
This is done for data from LEP and from various 
other e + e" colliders. Subsequently, the available 
results for £ * are compared to the LPHD + MLLA 
approach [1,2]. 

2 Experimental values of tl* for 
identified hadrons 

At LEP, analyses of explicitly identified hadrons 
have been performed for the ir° [6], the Kg [8,9], 
the r] [10] and the A [8,1 1,12,22]. At lower statistics 
data is also available for the E~ and fi~ baryons, 
as well as the S* and E* [8,11]. 

The values of correspond to low momenta, 
where the dependence on the primary quark flavours 
is expected to be small. To first approximation the 
distribution in £ p is gaussian, but a distorted gaus- 
sian fits the calculated spectrum more accurately 
[13,2]. The present statistics of identified hadrons 
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Figure 3: The maxima of the distributions in £ = — log x p for different hadrons as measured at LEP 
(table 1). the value of ^* is a nearly linear function of \ogQ . The exact dependence for A eff = 150 MeV 
leads to figure (b), where the naive expectation would be that Q is a linear function of the hadron mass. 
Figure (c) gives the dependence of Q on the hadron mass divided by the number of valence quarks. 



is not yet sufficient to extract the width and the 
two distortion parameters. However, since these 
other parameters depend more strongly on the high 
momentum data, they can be expected to depend 
stronger on the event flavour, which could bias the 
results. This is a reason not to use a distorted gaus- 

sian to determine f* . 

p 

In the following £* was determined by fitting the 
f p distribution in a limited range of f p around the 
maximum, to a gaussian distribution. A systematic 
error was estimated by changing the fitted range. 
The statistical and systematic errors were added in 
quadrature. The systematic error often dominated. 

The results are given in table 1 . Use was made 
of spectra published by CLEO [15], ARGUS [16], 
TPC [17], HRS [18], TASSO [19], JADE [20], 
TOPAZ [21], DELPHI [8,22], OPAL [5,9,11] and 
L3 [6,10]. Most of the recent publications reported 
values of £ * . Most experiments at lower centre of 
mass energies did not report the distribution of £ p or 
x p , but used the energy to define x E = 2E/E beam , 
or they divided by fi = v/c to obtain the 'scal- 
ing cross section'. Many different normalisations 
have been used, but fortunately this is unimpor- 
tant for a determination £ * . The experiments (AR- 
GUS and CLEO) near the T resonances give re- 
sults for the 'continuum', and for the resonances 
itself, where the contributions from the continuum 
are subtracted. The values in table 1 refer to data 
from the continuum above the T resonance. In the 
data from ARGUS for charged pions and protons, 



the contributions from decays of Kg and A were 
subtracted. 



3 Predictions about £ 



The LPHD + MLLA calculations [2] of the dis- 
tribution of f p depend on an effective QCD scale 
A eff ~ Aqcd and on a transverse momentum cutoff 
Qo in the evolution of the parton cascade. The value 
of £* is calculated to be a nearly linear function of 
logQo- 

The calculated dependence on the centre of mass 
energy is as follows [2] : 



Y 



2 

F(X) 



C/Y + 0(V 



-3/2- 



(1) 



where F(0) = and 1 



Y = log(£ beam /A eff ), A = log(Qo/A eff ) (2) 

and A eff is the 'effective' QCD scale, while only the 
cutoff scale Q depends on the hadron type. The 
constant C is calculated to be 



C 



4N r 



b ' 



(3) 



J This detail [23] is not completely clear in reference [2], 
where equation (1) is given as a function of Y — A, but in the 
approximation that A = 0. 
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where 

b = 11^/3-2^/3, (4) 

a = llN c /3 + 2n f /(3N*). (5) 

N c = 3 is the number of colours and n s = 3 is the 
active number of quark flavours in the fragmenta- 
tion process. For rif = 3 one finds C = 0.2915. 

It is important to note that F does not depend 
on Y, and that the first part of equation (1) is in- 
dependent of Q . This predicted behaviour can be 
checked by comparing spectra of different identi- 
fied particles. In the available momentum range 
this leads to a nearly linear dependence of on y. 

In ref. [2] various graphs show the results of 
numerical calculations for the dependence of vari- 
ous parameters on E be!im and Q ^ A eff . Figure 4 
of ref. [2] shows the dependence of on Qo for 
A eff = 150 MeV, and at a number of beam energies. 
From this it is possible to extract F after subtracting 
the Y dependent part of equation (1). For an inter- 
pretation of the experimental values of £ * it was 
useful to fit both this function F(X), and its inverse 
A(-F) to polynomials. The result for F(X) is: 

F(X) = -1.380-A + 0.086-A 2 + 0.044-A 3 ±0.06 

(6) 

while the inverse of F is described by: 

\(F) = -0.749- J F-0.054- J F 2 -0. 072-F 3 ±0. 06 

(V) 

The given accuracy is the maximum deviation from 
thedistributionplottedinref. [2], in the range —2 < 
F < and < A < 2. 

4 Dependence of ^* on the hadron 
mass 

The dependence of on the hadron mass as deter- 
mined from LEP data is shown in figure 3(a). Using 
equation (7) we can convert the values of for a 
given value of A eff = 150 MeV to values of Q . 
This result of this is shown in figure 3(b) 

Naively one could expect that the cutoff scale Q 
is larger for more massive hadrons. The most naive 
guess would be a linear dependence of the form 
<5o ~ A eff + m hadron . Figure 3(b) shows that this is 
not correct and that at least a separate treatment of 
mesons and baryons is necessary. A linear depen- 
dence of Q on e.g. the hadron mass per constituent 
quark already describes the data slightly better, as 
can be seen in figure 3(c). 




Figure 4: The experimental values of '£* versus \fs 
for various identified hadrons. The curves are from 
the combined fit of A e ff and F( hadron) to the pre- 
dicted dependence (equations (1) and (8)). From 
top to bottom, the curves correspond to pions ( it" 11 , 
ir°), all charged particles, kaons(K ± , K°), protons 
and A baryons. 

5 Dependence of on 

A more quantitative comparison is possible for the 
scaling behaviour of £*, as parametrised by the 
QCD scale A eff . The dependence on = 2_E bea m 
of the experimental values of £ * from table 1 is 
presented in the figure 4. 

The experimental values can be fitted easily to 
the dependence given by equation (1) with F = 
0, or to a linear dependence on Y. The results 
for different types of hadrons are given in table 
2. Both functions are acceptable and relation (1) 
gives approximately the correct slope, that should 
be equal for all hadrons. 

Before going on, one should note that in this case 
the quality of the fits can be considered as an internal 
check of the estimated errors. The 14 parameter 
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particle 


a-h 


b h 


X 2 /n.d.f. 
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7T± 


0.75 ±0.08 
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0.69 ±0.07 


0.715 ± 0.024 
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0.175 ± 0.004 


2.15 


14 


K± 


0.29 ± 0.11 


0.577± 0.035 


2.05 


0.589 ± 0.027 


2.24 


7 


K° 


0.46 ±0.12 


0.529 ± 0.035 


2.69 


0.612 ± 0.023 


5.91 


9 


K±,K° 


0.42 ±0.07 


0.540 ± 0.022 


2.14 


0.603 ± 0.018 


4.02 


16 


p 


0.33 ±0.15 


0.562 ± 0.046 


0.33 


0.621 ± 0.032 


1.07 


7 


A 


0.03 ±0.07 


0.610 ± 0.019 


0.34 


0.770 ± 0.025 


1.37 


7 


p,A 


0.12 ±0.07 


0.596 ± 0.018 


2.08 


0.725 ± 0.020 


2.08 


14 


total x 2 /n.d.f. 


1.35 


1.94 





Table 2: Result of a fit of the experimental values of in table 1 to a linear function = a h ± b h ■ 
log( E cms I GeV) for various identified hadrons. The scale Q' is the result of a second fit to the function 
= |y + V CY — C, where Y = log(E beam / Q' ) and C = 0.2915. The fits improve slightly when data 
at the T are excluded. To compare the overall quality of the fits to the two functions, a total x 2 /n.d.f. is 
given, based on the fits for ^ , ir°, K ± , K°, p, and A. 



fit to separate linear functions of log E for each 
hadron leads to a total x 2 /n.d.f. = 1.35. One could 
conclude that some of the (systematic) errors have 
been slightly underestimated. It can also be noted 
that the K° data give the largest contribution to x 2 - 

In the fit to equation (1) it is not easy to dis- 
entangle A eff from F(Qo), since they are strongly 
correlated. The simple way out was to constrain 
A eff = Q = Q' , so that F(X) = 0. The fitted 
values of Q' are not equal to the corresponding val- 
ues of Q shown in figure 3. The values of Q' are 
quite near to the particle mass for mesons, but for 
the baryons Q' is significantly smaller. 

Figure 5 shows the dependence of Q' on (a) the 
hadron mass and (b) the hadron mass per valence 
quark. Neither of these plots shows a clear linear 
dependence. It is possible that this is again due to 
F in equation (1). 

An alternative to setting F = is to perform a 
simultaneous fit of equation (1) to all the data from 
table 1 , where A eff is a universal scale, while F is 
allowed to have different values for hadrons, pions, 
kaons and the two baryons. 

Such a six parameter fit was performed by min- 
imising x 2 using the MINUIT [24] program, with 
zero as starting values for F h . There is a strong 
correlation in the fit of log A eff and F, the aver- 
age of F h . Therefore it is practical to fit F and 
differences F h — F. The fit has 44 degrees of free- 
dom and leads to x 2 /44 = 1.98, corresponding to 
a probability of 1.1 • 10~ 4 . This value of x 2 can 
also be compared to the values of x 2 /n.d.f in table 



2, where separate fits for every hadron type were 
done. This comparison shows that one should prob- 
ably not reject equation (1) on the basis of the low 
probability, because it cannot be excluded that the 
errors of £ * in table 1 have been underestimated by 
a factor of ^x 2 / n -d-f- = 1.41, on average. Figure 
4 shows that the fitted functions are very near to 
straight lines and that they follow the data points 
rather well. Addition of an 0(Y 3 l 2 ) term as in 
equation (1) does not improve the fit significantly. 

The resulting parameters and their external errors 
are given by: 

Aeff = 0.07 to 05 GeV 

F = -l.U+° f 

F w -F = 0.54 ±0.02 

^charged " F = 0.43 ± 0.02 (8) 

F K -F = -0.26 ±0.02 
F p -F_ = -0.28 ±0.03 
F A -F = -0.42 ±0.02 

The correlation coefficient of A eff and F is 100%, 
while the values of F h - F are constrained to have 
a zero sum. 

6 The theoretical meaning of A e ff 

The determined value of A eff from (8) is related to 
the (running) strong coupling constant by: 

^ = J_ = I (9 ) 

2tt bY Mog(£ beam /A eff )' 
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Figure 5: The values oj ? Q' from table from 2 versus 
(a) the hadron mass, and (b) the mass per valence 
quark. 

giving at the Z°: 

a s (LEP) = 0.108 ± 0.018 (10) 

This value is compatible with the more accurate 
determinations done e.g. at LEP. This would suggest 
that LPHD + MLLA can really give a quantitative 
description of the low x p data, based on perturbative 
QCD. 

However, theoretically the QCD scale A eff of 
MLLA is not well defined. This means that the 
value of A eff as determined in (8) can not be com- 
pared to known values of e.g. Aj-jg. A measurement 



of Aqcd hi a well defined subtraction scheme is 
equivalent to a measurement of a s (E besai ) at a well 
defined energy scale. To fix this energy scale one 
would need the (9(Y 3/,2 )-term in equation 1. 

Toseethis [23], onecan express!" in equation(l) 
in terms of the running coupling a s (E besm ). Using 
equation (9), the scale dependence of equation (1) 
translates to: 



d£*(-ET,eam, Qc 
d log E btma 



8N C 



' N c a s 
2tt 



o («r ) 

(11) 

The ^fa~ s term is the next-to-leading (MLLA) cor- 
rection to the leading term. The next-to-next-to- 
leading term 0{a s ) vanishes. An uncertainty of 
the energy scale would lead to e.g. a s (2E) instead 
of a s (E), inducing a correction in equation (11) 
of O (a 3 J 2 (E)). This shows that the O (a 3 J 2 (E)) 
needs to be known for the energy scale to be de- 
fined exactly. This is the reason for talking about 
the scale A eff , in stead of e.g. Aj-jg. 



7 Conclusion 

The data of LEP are giving a wealth of information 
about the fragmentation of quarks and gluons, as 
well as the mechanisms at play in the hadronisation. 
Although Monte Carlo models like JETSET and 
HERWIG can describe these data very accurately, it 
is questionable whether these models with all their 
free parameters lead to a better understanding of the 
physics behind the formation of hadrons in a jet, or 
only to a better parametrisation. 

The 'LPHD + MLLA approach [1-4] is an in- 
teresting attempt to understand some properties of 
multihadron production in jets in terms of perturba- 
tive QCD. This approach tries to take perturbative 
QCD as near as possible to the limits posed by 
confinement. It gives predictions for properties of 
momentum spectra that can be compared with ex- 
perimental data, but concentrating on observables 
that can be described without too many unknown 
free parameters. 

The experimental values of the parameter £* for 
various hadrons and at various centre of mass en- 
ergies have been determined from the published 
momentum spectra. Subsequently, an investigation 
was made of the dependence of £* on both the cen- 
tre of mass energy and the mass of the identified 
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hadrons, and this could be compared with the theo- 
retical predictions. 

The dependence of on the centre of mass en- 
ergy can be described adequately by the MLLA 
calculations. It is perhaps striking that the low x p 
data, that are so intimately related to confinement, 
can not only be described by perturbative QCD, but 
can even be used to extract a consistent value of the 
strong coupling constant. The dependence of £ * on 
the mass and flavour of the identified hadron raises 
some questions that have not yet been answered in 
the context of LPHD + MLLA and the description 
of fragmentation by truncated parton cascades. 
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